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ABSTRACT. 

[TJ shows, for a,kxk triangular matrix M = [rriij] with unique diagonal 
elements, how to obtain numbers Pij tS such that the (i, j) th element of M n 

3 

is given by n my = '^^Pij,s m ^ 1 ■ This paper generalizes that formula: 

s=i 

it removes the restriction of unique diagonal elements and shows how to 

num(i,j) mpy(r) , . 

obtain numbers Cij >r ,s such that „m, j = c ij',r,sf .. J m 7j,rJ 

r=l s=l ^ ' 

where num(i,j) is the number of unique diagonal elements between and 
including the i th and j th rows, {^ij^}""" 1 is the set of those unique 
elements, and mpy(r) is the multiplicity of mij >r on that same range. 
Like the Pi,j, s m [1]? the Cjj ; r,s are independent of the power to which the 
matrix is raised. The generalized formula works for any power of M, 
negative, zero or positive (positive only if the matrix is singular). 



Key words: Matrix, Triangular, Powers, Closed Form 



1. INTRODUCTION 

[1] presents a method of obtaining a simple closed form for the powers 
of a triangular matrix with unique diagonal elements, as follows: 

Definition 1. Let M = [rriij] be a k x k upper triangular matrix with 
unique diagonal elements. We define the power factors of M, Pij, s , re- 
cursively on the index j , as follows: 
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i-i 

Piis = — i<s<j<k, (1.1) 

Pi,j,s = 8<i, s> j, 

I>,,,,j m,,i - ^2 IJ '-! J i<j<k, (1.2) 

t=i 

Pj,j,j = m j,j- 

Theorem 1. If M = [rriij] is a non-singular upper triangular matrix 
with unique diagonal elements, and n m,ijis the (i,j) th element of M n , 
then 

3 

En— 1 
Pi,j,s m s ,s ' 

s=i 

for all integral values of n, negative, positive or zero. If M is singular, 
the equation holds if n > 1 and 0° is taken as 1). 



2. ALTERNATE DEFINITION FOR p i)jtS . 

Let M = [rriij] be an upper triangular matrix with unique diagonal 
elements. The product 

where i<a<b<c---<l<j, and s is i, j, or one of a,b,c,- ■ ■ I, is 
called a chain from i to j passing through s. The length of the chain is 
the number of elements in the product. The expression 

mi, a m a , b m b , c ■ ■ ■ rriij 

(m S!S - m aja )(m SjS - m feife )(m SjS - m CjC ) • • • (m S;S - m M )(m S;S - m jd ) 

where (m s>s — m s ^ s ) is taken as 1, is called an adjusted chain from i to j 
passing through s. 
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Definition 2. If i < s < j, Pij, s is the sum of all adjusted chains from 
i to j passing through s. If s < i or s > j, Pi j S = 0. 

Following are a few illustrative examples which help clarify the defi- 



mi,imi,3 m 1 , 1 m 1 . 2 m 2 ,3 

(mi,i-m3,3) (mi,i— m2,2)(mi,i— m,3,3), 

"ti i 2"i2,3 j rai.irai, 2^2,3 

m.2,2— "13,3 (m2,2— mi,i)(rn.2,2— m.3,3), 

I mi,imi,3 _|_ ?»i,2"i2,3 j rn ltl m li2 m 2 ,3 

1,3 "13,3—011,1 m3 : 3-rrt2,2 ("13,3— mi,i)( m 3,3-rn2,2)- 

Theorem 2. Definition 2 is equivalent to Definition 1. 

Proof. From Definition 2, each term of the summand in (1.1) is of the 
form 

mi t ama tb m btC ■ ■ ■ m^ t m td 

(ra s , s - m ata )(m StS - m b)b )(m S)S - m c>c ) • • • (m S)S - m M )(m SjS - m t ,t) ' 

where (m SjS — m SjS ) is taken as 1. The sum of all such terms, from t = s 
to t = j — 1, clearly includes all of the chains from i to j passing through 
s. That sum would be the sum of all the adjusted chains from i to j 
passing through s, except that the difference (m SjS — rrij ^ is missing from 
each denominator. Hence the division by that difference in (1.1), and 
thus (1.1) is satisfied. 

j 

Next, we need to show that (1.2) is satisfied. The sum Pi,j,t consists 

t=i 

of all adjusted chains from i to j of the form 

'^i,a\ n ^a\,a 2 P^ J a 2 ,az ' ' '^a r ,j 

(m S!S - m aiiai )(m 3ia - m a ^ a2 ){m s ^ s - m a3ia3 ) ■ • • (m s>s - m ar ^ ar ){m SiS - m jd ) ' 

with values r = 1, 2, 3, • ■ ■ , j — i and s = a\, 02, 03, • ■ • ,j . The only 
term with r = 1 comes from pijj and is equal to m^j (see the illustrative 



nit ion: 

Pi, 1,1 = 

Pl.3.1 = 
Pl.3.2 = 
Pl,3,3 = 
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example for p^s,3 in Definition 2). The sum of all the adjusted chains of 
length 2 with the same numerator is 

m^.m,,, a A 1 ) = 0. 

The sum of all adjusted chains of length 3 with the same numerator 

is 

m i,a 1 m a 1 ,a2 m a 2 ,a 3 ( ( ma ^ ^ -m a2>a2 )(m al ,^-"^3,03) {m a2 , a2 -m alial )(m a2> a 2 -m a3>a3 ) ^ 
7 k t) = 0. 

In general, the multiplier of mi ax m ai ^ 2 m a2AZ ■ ■ ■ m ar j is seen [2] to be 
the (r — l) st divided difference of the polynomial f(x) = 1, and hence is 
if r > 2. 
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Therefore, Pi,j,t — m i,j > an d (1-2) is satisfied. And since Pj, h3 



t=i 

m it j, Definition 2 is equivalent to Definition 1. 



□ 



3. NON-UNIQUE DIAGONAL ELEMENTS 

Theorem 3. Let M = [rriij] be an upper k x k triangular matrix with 
non-unique diagonal elements. Let num(i,j) be the number of unique di- 
agonal elements between and including the i th and j th rows, {jrij^r}™"" 1 
be the set of those unique elements, and mpy{r) be the multiplicity of 

num(x,y) 

rriij :r on that same range ( mpy{r) = k). 

t=i 

Then 

num(i,j) mpy(r) , _^ 



r=l s=l 
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where the Cij, r ,s are independent of the power to which the matrix is raised. 
This generalized formula works for any power of M , negative, zero or 
positive (positive only if the matrix is singular). 



Proof. The proof is much easier to present, and much easier to follow, 
for a specific case. The generalization in Theorem 3 will be evident, 
including the method of determining the numerical values of c%,j,r,8- 



Let 



M 



( 3 2 3 5 4 

5 2 4 3 

3 2 6 

5 5 

7 

\ 



2\ 
1 
4 
1 

2 

3 / 



In terms of Theorem 3, 



mi, 6| i = 3 mpy{\) = 3, 
"7,1,6,2 = 5 mpy(2) = 2, 
mi,6,3 = 7 mpy(3) = 1. 

and num(l, 6) = 3. In order to make use of Theorem 1, we first alter the 
matrix M so that it will have unique diagonal elements, as follows: 



M 



( 3 2 

5 







y 



3 


5 


4 


2 


2 


1 


3 


1 


— e 


2 


6 


1 





5-/ 


5 


1 








7 


2 











3 - 



9 J 



At an appropriate stage, we will let e, /, and g become zero to obtain 
results for the original matrix M. 
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We have from Theorem 1, for example, the element of M n at row 1 
and column 6 is given by 



n mifi = Pi.6,13"- 1 + Pi, 6| 25 n " 1 + pi, 6 ,3(3 - e)"- 1 + (3.1) 

Pi,6,4(5 - fT- 1 + PififiT*- 1 + p 1)6 , 6 (3 - g)"- 1 . 

Expanding the binomials, and collecting only those terms involving 
powers of 5, we obtain 



n-l , _ i \ 

(Px^ + Pw^+Pw^f . ){-m M . (3.2) 

i=l ^ % ' 

We know from Definition 2 that pi,6,4 is the sum of all adjusted chains 
from 1 to 6 passing through 4. Combine all of those adjusted chains 
into a single fraction with the denominator equal to the product of the 
differences between m 4j4 and each of the other diagonal elements, namely 
(2-/)(-/)(2 + e-/)(-2-/)(2 + £-/). We see from this denominator 
that Pi,6,4( — fY = if i > 1 and / = 0. Therefore, the only terms of 
(3.2) involving powers of 5 are 

(Pl,6,2 + Pl,6 i 4)5"" 1 + Pl, M (-/) ^ ~ ^ 5"" 2 . 

We determine the coefficients of the powers of 5, Ci i6i2 ,i and Ci j6j 2,2, 
using Mathematica as follows: 

ci, 6>2)1 Input: (p[l,6,2]+p[l,6,4]//Together)/.{e- > 0, /- > 0,g- > 

0} 

Output: -f 

c 1A2 ,2 Input:((-/)p[l,6,4]//Together)/.{e- > 0, /- > 0,0- > 0} 
Output: —60 
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Similarly, collecting only those terms involving powers of 3, we obtain 



(Pl,6,l+Pl,6,3+Pl,6,6)3 n " 1 + ^(pi, 6 ,3(-e) i +Pl,6,6(-^) 1 ) . 3 




(3.3) 



Combine all of the adjusted chains in (pi,6,3(— e) 1 + Pi$,q{— g)Y into a 
single fraction, with the denominator equal to a x b, where a is the 
product of the differences between 771.3,3 and each of the other diagonal 
elements, and b is the product of the differences between m^^ and each 
of the other diagonal elements. The only elements in that denominator 
which become zero when e, /, and g become zero are (— e)(g — e)(e — 
9){-g) = eg(g-e) 2 . If {pi e) 1 +pi fifi (-g) % ) were equal to a constant 
A 7^ 0, the numerator of that single fraction would have to contain a 
term eg(g — e) 2 A. If i > 2, that is not possible since each term of the 
numerator must contain a power of e greater than 2 or a power of g 
greater than 2, but eg(g — e) 2 A contains the term (—2)e 2 g 2 A. Hence 
(pifi,3(— e) 1 +Pi,6,e)(— g) 1 = if i > 2 and e, /, and g are zero. Therefore 
the only terms of (3.3) involving powers of 3 are 



Again, we determine the coefficients of powers of 3, c 1)611 , Ci i6j i )2 , Ci )6) i j3 , 
using Mathematica as follows: 



ci j6jl ,i Input: (p[l, 6, 1] + 6, 3] + p[l, 6, 6] //Together) /.{e- > 0, /- > 0, g- > 0} 
Output: 

Cl , 6)1>2 Input: ((-e)p[l,6,3] + {-g)p[l, 6, 6]//Together)/.{e- > 0,/- > 0,^- > 0} 




Output: I 
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ci )6il ,3 Input: ((-e) 2 p[l,6,3] + {-g) 2 p[l, 6, 6] //Together) /.{e- > 0, /- > Q,g- > 0} 
Output: y 

The only term in (3.1) involving powers of 7 is Pifi^7 n ~ l . We de- 
termine the coefficient of 7 n_1 , 01,6,3,1, using Mathematica, as follows: 

ci, 6 ,3,i Input: (p[l,6,5]//Together)/.{e- > 0,/- > 0,g- > 0} 
Output: 

Combining the above results we obtain the closed form expression for 



It is easy to see from the pattern in (3.1) that the structure of Theo- 
rem 3 is correct in the general case. To prove Theorem 3 we need only 
establish the cutoff points for the expansion of the binomials in (3.1). 
We did that explicitly for elements of multiplicity 2 and 3. We do so now 
for elements of higher multiplicities. 

The argument immediately following (3.3) proved that for a diagonal 
element of multiplicity 3, the expression (pi,6,3(— e) 4 + pifi,e(—g) 1 ) was 






m l,6,l + C l,6,l 
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zero if i > 2 and e and g were zero. It did that by showing that the 
numerator would have to contain a term eg(g — e) 2 A(A ^ 0), which was 
then shown to be not possible. 

Suppose now that the element m 5j5 in the unaltered matrix M were 
3 instead of 7, and was replaced by 3 — h in the altered matrix. That 
same argument, but now for a diagonal element of multiplicity 4, would 
change as follows: 

(1) (pi,6,3(-ey+p lfifi (-gY) would be replaced by {pi^{-e) l +p lfifi {-h) l + 

PififiirgY), 

(2) eg(g — e) 2 would be replaced by egh(e — g) 2 (e — h) 2 (g — h) 2 , 

(3) If (pi,6,3(~e)* +Pi,6,5( — h) % + Pififi(— g) 1 ) were equal to a constant 
A ^ 0, the numerator would have to contain a term egh(e — g) 2 (e — 
h) 2 (g — h) 2 A. But if % > 3, that would not be possible because each term 
of the numerator would have to contain a power greater that 3 for at 
least one of e,g, or h, but egh(e — g) 2 (e — h) 2 (g — h) 2 A contains the term 
-8e 3 g 3 h 3 A. 

The generalization to higher multiplicites is straightforward. 

□ 



4. ILLUSTRATIONS 

(1) We use Theorem 3 to obtain a closed form for n m,2^ for the ma- 
trix M in Section 3, noting that 7712,4,1 = 5, vnpy{\) = 2, 772-2.4,2 = 3, 
mpy(2) = 1, and num(2,4) = 2. The coefficients C2^, r .s, are determined 
using Mathematica as follows: 

c 2 ,4,i,i Input: (p[2,4,2] + p[2, 4, 4]//Together)/.{e- > 0, /- > 0} 
Output: 1 
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c 2AX2 Input:((-/)p[2,4,4]//Together)/.{e- > 0, /- > 0} 
Output: 30 

02,4,2,1 Input:(p[2,4,3]//Together)/.{e- > 0, /- > 0} 
Output: 3 

Theorem 3 gives us 

n m 2i4 = 1 • 5 n ~ x + 30 ~ l \ 5 n ~ 2 + 3 ■ 3"" 1 . 

(2) Let M and the altered M be the matrices 

/ 5 2 1 3 \ / 5 2 1 3 \ 

5 4 2 5-e 4 2 

5 1 5-f 1 

\ 5 / \ 5-g J 

We use Theorem 3 to obtain a closed form for n mi,4 , noting that 
"T-i,4,i = 5, mpy(l) = 4, and ratm(l,4) = I. The coefficients 
determined using Mathematica as follows: 

ciah Input: (p[l,4,l] + p[l,4,2] +p[l,4,3] +p[l,4,4]//Together)/.{e- > 0, /- > 0},^- > 0} 
Output: 3 

ci,4,i, 2 Input:((-e)p[l,4,2] + {-f)p[l, 4, 3] +p{-g)p[l, 4, 4]//Together)/.{e- > 0, /- > 0,0- > 0} 
Output: 20 

ci, 4 ,i,3 Input:((-e) 2 p[l,4,2] + (-/) 2 p[l, 4, 3] + p(-^) 2 p[l, 4, 4]//Together)/.{e- > 0, /- > g- > 0} 
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Output: 33 



ci.4,1,4 Input:((-e) 3 p[l,4,2] + (-/) 3 p[l, 4, 3] + p(-^) 3 p[l, 4, 4]//Together)/.{e- > 0, /- > g- > 0} 
Output: 40 

Theorem 3 gives us 
n m M = 3 • 5"" 1 + 20 ( H ~ l \ 5"~ 2 + 33 ( U ~ X ) 5^ 3 + 40 ( U ~ *) 5"" 4 . 



Note: In the above illustrations, three different variables (e,f and g) 
were chosen for clarity of presentation. In actual practice, replacing (e,/ 
and g) with variables such as e, e 2 , and e 3 would be more efficient since 
Mathematica would need to record only one variable instead of three, 
and only one variable need be set to zero. 
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